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Abstract. 

It is well known that one can map certain properties of random matrices, fermionic 
gases, and zeros of the Riemann zeta function to a unique point process on the real 
line M. Here we analytically provide exact generalizations of such a point process in 
d-dimensional Euclidean space R'* for any d, which are special cases of determinantal 
processes. In particular, we obtain the n-particle correlation functions for any n, 
which completely specify the point processes in M.'^. We also demonstrate that spin- 
polarized fermionic systems in R'^ have these same n-particle correlation functions in 
each dimension. The point processes for any d are shown to be hyperuniform, i.e., 
infinite wavelength density fluctuations vanish, and the structure factor (or power 
spectrum) S{k) has a nonanalytic behavior at the origin given by S{k) ^ |fc| (fc — + 0). 
The latter result implies that the pair correlation function 52 (j") tends to unity for 
large pair distances with a decay rate that is controlled by the power law l/r''+^, 
which is a well-known property of bosonic ground states and more recently has been 
shown to characterize maximally random jammed sphere packings. We graphically 
display one- and two-dimensional realizations of the point processes in order to vividly 
reveal their "repulsive" nature. Indeed, we show that the point processes can be 
characterized by an effective "hard-core" diameter that grows like the square root of d. 
The nearest-neighbor distribution functions for these point processes are also evaluated 
and rigorously bounded. Among other results, this analysis reveals that the probability 
of finding a large spherical cavity of radius r in dimension d behaves like a Poisson 
point process but in dimension d+1, i.e., this probability is given by exp[— K(c?)r''+"'^] 
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for large r and finite d, where K(d) is a positive rf-dependent constant. We also show 
that as d increases, the point process behaves effectively like a sphere packing with 
a coverage fraction of space that is no denser than 1/2''. This coverage fraction has 
a special significance in the study of sphere packings in high-dimensional Euclidean 
spaces. 
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1. Introduction 

It is well known that there is a remarkable connection between the statistical properties 
of certain random Hermitian matrices, the zeros of the Riemann zeta function, and 
fermionic gases [U [Sj El HI IH [6]. Underlying each of these objects are certain one- 
dimensional point processes (defined more precisely in Section 2) whose statistical 
properties (under certain limits) are believed to be identical. The purpose of this 
paper is to provide generalizations of this unique point process to point processes in 
(i-dimensional Euclidean space M.'^ for any d, and to characterize their spatial statistics 
exactly. Since the characterization of a point process can be viewed as the study of a 
system of interacting "point" particles, exact descriptions of nontrivial point processes 
in arbitrary space dimension, which are hard to come by [5, are of great value in the 
field of statistical mechanics. 

There are three prominent theories of random Hermitian matrices, which model the 
Hamiltonians of certain random dynamical systems; see the excellent book by Mehta 
[3]. If the dynamical system is symmetric under time reversal, then the relevant theory 
for integral spin is that of the Gaussian orthogonal ensemble (GOE) or the Gaussian 
symplectic ensemble (GSE) for half-integer spin. On the other hand, the Gaussian 
unitary ensemble (GUE) models random Hamiltonians without time reversal symmetry, 
which is relevant to certain properties of the Riemann zeta function. Although there 
are distinct one-dimensional point processes associated with each of these ensembles, 
our interest here will be in the GUE because of its relationship to the Riemann zeta 
function. 

The GUE of degree consists of the set of all x Hermitian matrices together 
with a Haar measure. This is the unique probability measure on the set of N x N 
Hermitian matrices that is invariant under conjugation by unitary matrices such that 
the individual matrix entries are independent random variables. Dyson [1] showed that 
the eigenvalue distrubutions of the GUE are closely related to those of the "circular 
unitary ensemble" (CUE), which he exactly mapped into a problem of point particles 
on a unit circle interacting with a two-dimensional Coulombic force law at a particular 
temperature. This point process on the unit circle or, equivalently, on the real line M in 
the large- A^ limit (when suitably normalized) has a pair correlation function (defined in 
Section 2) in M at number density p = 1 given by 

, , ^ sin^(7rr) 

Equation ([T]) also applies for the GUE in the limit A^ +oo such that the mean gap 
distance between eigenvalues at the origin is normalized; this limit has the effect of 
magnifying the bulk of the eigenvalue density on M such that ([T]) is well-defined. We 
see that this point process is always negatively correlated, i.e., g2{r) never exceeds unity 
(see Fig. [1]) and is "repulsive" in the sense that g2{r) ^ as r tends to zero. More 
generally, Dyson P| proved that the n-particle correlation function (defined in Section 
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2) is given by the following determinant: 



«J / ij=l,...,n 
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Figure 1. Left panel: The pair correlation function g2{r) versus distance r for the 
eigenvalues of the GUE/CUE in the large- limit, which is conjectured to be the same 
as the one characterizing the nontrivial zeros of the Riemann zeta function. Spin- 
polarized fermions in their ground state in M have the same pair correlation function. 
Right panel: The corresponding spectral counterpart, called the structure factor S{k) 
[cf. ([7])], as a function of wavenumber k. 

Montgomery [3] conjectured that the pair correlation function of nontrivial zeros 
of the Riemann zeta function (which, according to the Riemann hypothesis, lie along 
the critical line 1/2 + it with t G M) is given exactly by the GUE/CUE function ([1]) 
in the asymptotic limit (high on the critical line) when appropriately normalized. This 
remarkable correspondence was further established by Odlyzko [5], who numerically 
verified the Riemann hypothesis for the first 10^^ nontrivial zeros of the zeta function as 
well as at much larger heights and confirmed that the pair correlation function agrees 
with ([1]). Rudnick and Sarnak [6] proved that, under the Riemann hypothesis, the 
nontrivial zeros have n-particle densities for any n given by ([2]). The reader is referred 
to the excellent review article by Katz and Sarnak [7], which discusses the connection 
between the zeros of zeta functions and classical symmetric groups, of which the three 
canonical random-matrix ensembles are but special cases. 

For spin-polarized free fermions in M (fermion gas) at number density p = 1, it is 
known that the pair correlation function in the ground state (i.e., completely filling the 
Fermi "sphere" ) is given by 

where kp is the Fermi radius, which is the one- dimensional analog of the well-known 
three-dimensional result [1]. Therefore, we see that when kp = vr, we obtain the 
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GUE/CUE pair correlation function ([T]). The repulsive nature of the point process 
in this context arises physically from the Pauli exclusion principle. 

Very little is known about how to generalize one-dimensional point processes 
associated with random matrices and number-theoretic functions to higher dimensions. 
In this paper, we analytically obtain exact generalizations of the aforementioned one- 
dimensional point process in dimensional Euclidean space for any d. These 
processes are special cases of determinantal processes. In particular, we obtain the 
n-particle correlation functions for any n, which completely specify the point processes 
in M."^. We also show that spin-polarized fermionic systems in have these same n- 
particle correlation functions in each dimension. We show that the point processes for 
any d are hyperuniform, i.e., infinite wavelength density fluctuations vanish, and the 
structure factor (or power spectrum) S{k) has a nonanalytic behavior at the origin 
given by S{k) ~ \k\ {k ^ 0). The latter result implies that the pair correlation function 
g2{r) tends to unity for large pair distances with a decay rate that is controlled by 
the power law In three dimensions, such a dominant power- law decay of g2{r) 

is a well-known property of bosonic systems in their ground states [TOl [I2] and, more 
recently, has been shown to characterize maximally random jammed sphere packings 

m- 

Realizations of the point processes are displayed in one and two dimensions, 
using a simulation technique described by us elsewhere [H], which vividly reveal their 
"repulsive" nature. In fact, we show the point processes can be characterized by 
an effective "hard-core" diameter that grows like the square root of d. The nearest- 
neighbor distribution functions for these point processes are also studied by evaluating 
and rigorously bounding them. Among other results, this analysis reveals that the 
probability of finding a large spherical cavity of radius r in dimension d behaves 
like a Poisson point process but in dimension d + 1, i.e., this probability is given by 
exp[— K,{d)r'^^^] for large r, where K{d) is a positive d-dependent constant. We also show 
that as d increases, the point process behaves effectively like a sphere packing with a 
coverage fraction of space that is no denser than 1/2'^. As we will see, this value of 
the coverage fraction has a special significance in the study of sphere packings in high 
dimensions. 

In Section 2, we present background and definitions concerning point processes 
that are particularly germane to the ensuing analysis. Section 3 discusses general 
determinantal point processes. In Section 4, we obtain the determinantal point processes 
in R'^ that are generalizations of the aforementioned one-dimensional point process 
associated with fermions, random matrices, and the Riemann zeta function. We call the 
most general of these point processes "Fermi-shells" point processes. The asymptotic 
properties of various pair statistics are investigated. We show that spin-polarized 
fermionic systems in M'^ have the same n-particle correlation functions in each dimension. 
Section 5 analyzes various nearest-neighbor functions for the special case of the "Fermi- 
sphere" point processes in R"^. We present concluding comments in Section 6, including 
remarks about possible connections of our point processes to random matrix and number 
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theory. 

2. Background on Point Processes 
2.1. Definitions 

A stochastic point process in M'^ is defined as a mapping from a probabihty space 
to configurations of points xi,X2,X3... in ci-dimensional Euchdean space M.'^. More 
precisely, let X denote the set of configurations such that each configuration x E X is 
a subset of M*^ that satisfies two regularity conditions: (i) there are no multiple points 
(xj 7^ Xj a i ^ j) and (ii) each bounded subset of M*^ must contain only a finite number 
of points of X. We denote by N{B) the number of points within x (1 B, B E B, where 
B is the ring of bounded Borel sets in W^. Thus, we always have N{B) < oo for B E B, 
but the possibility A^(R'^) = oo is not excluded. We note that there exists a minimal 
cr-algebra U of subsets of X that renders all of the functions N{B) measurable. Let 
JF, V) be a probabihty space. Any measurable map x : Q ^ X is called a stochastic 
point process ^B, [ISj . Henceforth, we will simply call this map a point process. Note 
that this random setting is quite general. It incorporates cases in which the location of 
the points are deterministically known, such as a lattice. 

A point process is completely statistically characterized by specifying the 
countably infinite set of n-particle probability density functions Pn(ri, r2, . . . , r„) (n = 
1,2,3...) [16]. The distribution-valued function p„(ri, r2, . . . , r„) has a probabilistic 
interpretation: apart from trivial constants, it is the probability density function 
associated with finding n different points at positions ri,...,r„ and hence has the 
nonnegativity property 

p„(ri,r2,...,r„) >0 Vr, G M'^ = 1,2, . . .n). (4) 
The point process is statistically homogeneous or translationally invariant if for every 
constant vector y in M'^, Pn(ri, i'2, • • • , fn) = Pn(ri + y, • • • , + y), which implies that 

p„(ri, r2, . . . , r„) = p"5'„(ri2, . . . , ri„), (5) 

where p is the number density (number of points per unit volume in the infinite- volume 
limit) and 5'„(ri2, . . . ,ri„) is the n-particle correlation function, which depends on the 
relative positions ri2, ri3, . . ., where r^j = r^ — rj and we have chosen the origin to be at 
ri. We call 5'2(r) = 5'2(~r) the pair correlation function. 

For translationally invariant point processes without long-range order, 5'n(i'i2, • • • , fin) 
1 when the points (or "particles") are mutually far from one another, i.e., as \rij\ oo 
(1 < i < j < oo), Pn(ri, ^2, . . . , r„) — > p". Thus, the deviation of gn from unity provides 
a measure of the degree of spatial correlation between the particles with unity corre- 
sponding to no spatial correlation. Note that for a translationally invariant Poisson 
point process, gn is unity for all values of its argument. 

It is useful to introduce the total correlation function h{r) of a translationally 
invariant point process, which is related to the pair correlation function via 

h{r)=g2{r)-l (6) 
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and decays to zero for large |r| in the absence of long-range order. Note that h{r) = 
for all r for a translationally invariant Poisson point process. An important nonnegative 
spectral function S(k), called the structure factor (or power spectrum), is defined as 
follows: 

S{k) = l + phik), (7) 

where /i(k) is the Fourier transform of h{r). For a translationally and rotationally 
invariant point process, it is useful to consider the cumulative coordination number 
Z{r), defined to be the expected number of points found in a sphere of radius r centered 
at an arbitrary point of the point process, which is related to the pair correlation function 
as follows: 

Z{r)=psi{l) [ x''''g2{x)dx, (8) 



where 







d-1 



is the surface area of a d-dimensional sphere of radius r. It is clear that since g2{f') is a 
nonnegative function, Z{r) is a monotonically increasing function of r. 
The Fourier transform of an function / : R'^ ^ M is defined by 

/»= / /(r)exp[-i(k-r)]dr, (10) 

where denotes the space of absolutely integrable functions on M^. If / : M'^ ^ M is 
a radial function, i.e., / depends only on the modulus r = |r| of the vector r, then its 
Fourier transform is given by 



oo 



/(*) = (2.)t r'-^nn^^m^dr, (H) 

where k is the modulus of the wave vector k and Ju{x) is the Bessel function of order 
u. The inverse transform of f{k) is given by 

rk--^m'-ji:^dk. (12) 

(27r)2 Jo {kry'^^^> ^ 

2.2. Number Variance and Hyperuniformity 

We denote by o''^{A) the variance in the number of points N{A) contained within a 
window A C M'^. The number variance cr^iA) for a specific choice of A is necessarily 
a positive number and is generally related to the total correlation function h{r) for a 
translationally invariant point process [TS]. In the special case of a spherical window of 
radius R in M.'^, it is explicitly given by 



a\R)=pvi{R) 



I + p h(r)a(r; R) dr 



(13) 



where (J^{R) is the number variance for a spherical window of radius R in and a(r; R) 
is the volume common to two spherical windows of radius R whose centers are separated 
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by a distance r divided by vi{K), the volume of a spherical window of radius R, given 
explicitly by 

We will call a(r; R) the scaled intersection volume. 

For large R, it has been proved that (t'^{R) cannot grow more slowly than 'yR'^~^, 
where 7 is a positive constant [Hj. We note that point processes (translationally 
invariant or not) for which cr^(i?) grows more slowly than the window volume (i.e., 
as R'^) for large R are examples of hyperuniform (or superhomogeneous) point patterns 
[IHlIin]. This classification includes all periodic point processes [TH|, certain aperiodic 
point processes [HI [18], one-component plasmas [HI [18], point processes associated 
with a wide class of tilings of space [SUl [21], and certain disordered sphere packings 
[m [221 [23], [21] . Hyper uniformity implies that the structure factor S(k) has the following 
small k behavior: 

lim S(k) = 0. (15) 

The scaled intersection volume a{r; R) appearing in ( |T3|) and its associated Fourier 
transform will play a prominent role in this paper. The former quantity is defined by 
a convolution integral involving the indicator function w for a d-dimensional spherical 
"window" of radius R centered at position xq [E], i.e.. 



^i;(|x-xo|;i?) = 0(i?-|x-xo|), (16) 
where B(x) is the Heaviside step function 

Specifically, the scaled intersection volume is given by 

a(r; R) = ] / w(ri - Xq; R)w(r2 - Xq; R)dxo. (18) 

^^l(^) jRd 

The scaled intersection volume has the support [0, 2R], the range [0, 1], and the following 
alternative integral representation [T6]: 

/.cos-i[r/{2R)] 

a{r;R) = c{d) / sin'^(0) d^, (19) 







where c{d) is the dimensional constant given by 

2Vil + d/2) , , 



Torquato and Stillinger [16] found the following series representation of the scaled 
intersection volume a{r\R) for r < 2R and for any d: 

air; R) = I - c{d)x + c{d) y^_i)n {d-m-^)---{d-2n + 3) 

' ^ (2n-l)[2-4-6---(2n-2)] ' ^ ' 
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where x = r/{2R). For even dimensions, relation fl21l) is an infinite series, but for odd 
dimensions, tlie series truncates sucli tliat a(r; R) is a univariate polynomial of degree 
d. In even dimensions, the scaled intersection volume involves transcendental functions 
(e.g., for d = 2, a{r;R) = 27r-i[cos-i(27j) - - j^Y^^] for < r < 2R). Figure [2] 
provides plots of a(r; R) as a function of r for the first five space dimensions. For any 
dimension, a(r; R) is a monotonically decreasing function of r. At a fixed value of r in 
the interval (0, 2R), a(r; R) is a monotonically decreasing function of the dimension d. 

The Fourier transform of the scaled intersection volume function (HM . which is 
given by 



^ik;R) = \jj^ r^/V(./2M(A;r)dr 



where 

infh- m — 

271 \ '^1'^ 

—J R'J,/2{kR), (23) 

is the Fourier transform of the window indicator function f|T6|) . Therefore, the Fourier 
transform of a{r; R) is the following nonnegative function of k: 

a{k- R) = 2V/^Til + d/2) (^^^^) ' . (24) 




Figure 2. Left panel: The scaled intersection volume a{r; R) for spherical windows 
of radius i? as a function of r for the first five space dimensions. The uppermost curve 
is for d — 1 and lowermost curve is for d = 5. Right panel: Corresponding Fourier 
transforms for the case R = 1. The uppermost curve is for d = 5 and lowermost curve 
is for d ~ 1. 

It has been shown that finding the point process that minimizes the number variance 
(J^{R) is equivalent to finding the ground state of a certain repulsive pair potential with 
compact support [18]. This problem is directly related to an outstanding problem in 
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number theory involving generalized zeta functions and lattices Understanding 
such ground states can be facilitated by utilizing "duality" relations that link ground 
states in real space to those in Fourier space [26] . 

3. Determinantal Point Processes 

We will be able to obtain the appropriate d-dimensional generalizations of point 
processes corresponding to the eigenvalues of the GUE, the zeros of the Riemann 
zeta function, or the one-dimensional fermionic gas by appealing to the notion of a 
determinantal point process in M'^ [27]. Determinantal point processes were introduced 
by Macchi [27], who originally called them fermion point processes. Soshnikov [28] 
presented a review of this subject and discussed applications to random matrix 
theory, statistical mechanics, quantum mechanics, and representation theory. It is also 
noteworthy that examples of determinantal point processes arise in self-avoiding random 
walks [30] and uniform spanning trees [3T] . 

Without loss of generality, the number density is set to unity (p = 1) in the ensuing 
discussion. Let H{r) = H{—r) be a translationally invariant Hermitian-symmetric 
kernel of an integral operator Ti. A translationally invariant determinantal point process 
in exists if the the n-particle density functions for n > 1 are given by the following 
determinants: 

Pn(ri2, ri3, • • • , rin) = det[if(rij)]ij=i,...,„, (25) 

where H[0) = 1. By virtue of the nonnegativity of the p„ in the pointwise sense [cf. dl])] 
and ( l25l) . it follows that Ti must have nonnegative minors and Ti. must be a nonnegative 
operator, which implies that H{r) is positive semidefinite. The latter implies that the 
Fourier transform H{k) of the kernel H{r) is nonnegative, and this property together 
with the condition H{0) = 1 = J^^ if(k)(ik implies that H{k.) < 1, i.e., 

< //(k) < 1 for all k. (26) 

It follows that any positive semidefinite Hermitian-symmetric kernel H{r) = H{—r) 
whose Fourier transform satisfies the inequalities in ( l26l) describes a determinantal point 
process with a pair correlation function given by 

g2iv) = l-\Hir)\', (27) 

such that 

< ^2(r) < 1 and c/2(0) = 0, (28) 

and a n-particle density given by ( l25l) . We see that the total correlation function for a 
determinantal point process is given by 

h{r) = -\H{r)\\ (29) 

The fact that the n-particle density functions can be written in terms of the 
determinant specified by fl25l) leads to bounds on p„. For example, it trivially follows 



Point processes in arbitrary dimension 



11 



that 

Pn(ri2,ri3, • • • ,rin) < 1- (30) 
A less obvious but stronger upper bound is as follows: 

Pn(ri2, ri3, . . . , ri„) < P2(ri2)p2(ri3) • • ■ P2(ri„). (31) 
We remark also that: 

p„(ri2, . . . ,ri„) < pm(ri2, . . . ,rim) Vm < n. (32) 

Each of these inequalities is a consequence of the determinantal form for p„ and the 
characteristics of H; in particular, ( l32i) follows directly from Fischer's inequality [32] 
and an appropriate partition of the matrix representation for the operator 7i. Equation 
(pTl) is a result of the more general Hadamard-Fischer inequalities [32] and the relation 
if (0) = 1. The positive semidefinite character of H is essential for these inequalities to 
hold. 

A unique property of determinantal point processes is that each of the n-particle 
correlation functions Qn can be expressed completely in terms of the pair correlation 
function g2. Namely, at unit density we may write: 

g^ir) = 1 - [Hir)f (33) 

^ H{r) = ±v/l-(72(r). (34) 

Therefore: 

^„(ri2, . . . , ri„) = det [H{rij)]^^.^.^^ 
= det 



±^l-g2{r,,) . (35) 

J l<i<j<n 

The right side of ( l34l) is well-defined for all r G since < g2{r) < 1 Vr. We note 
that in general H{r) may be either positive or negative for a given value of r as in ([2]); 
therefore, the sign of the square root in (IMI) must be chosen appropriately. Our ability to 
express each gn in terms of the pair correlation function g2 is a reflection of the fact that 
the 77,-particle correlation functions depend on a common kernel H; such a reformulation 
is generally not possible for an arbitrary point process. Thus, one can infer the behavior 
of the n-particle correlation functions for n > 3 solely from a knowledge of the behavior 
of g2. 

A trivial example of a determinantal point process is the case in which H{0) = 1 
and H{r) = for r ^ 0. The resulting pair correlation function is given by g2{0) = 
and g2{r) = 1 Vr 7^ 0, and this function belongs to the same equivalence class as the 
pair correlation function of the Poisson point process with respect to Lebesgue measure 
on the nonnegative reals M"*" (i.e., the functions differ only on a set of measure zero). 
Note that Costin and Lebowitz [35] have considered the conditions under which a pair 
correlation function of the form g2{r) = 1 — exp(— Ar) is a determinantal point process 
in M^. 

We note in passing that the number of points in a determinantal point process 
that falls in a compact set A C M'^ has the same distribution as a sum of independent 
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Bernoulli(A^) random variables, where are the eigenvalues of the operator Ti, 
Moreover, Hough et al. [29] presented an algorithm to generate determinantal point 
process in W^, which we apply elsewhere [1] 



92{r) = 1 - ^7TZ^^■ (37) 



4. New Determinantal Point Processes in and Their Connection to 
Fermionic Gases 

"Fermi-Sphere" Point Processes in 

Here we obtain the appropriate generalization of ([1]) that corresponds to a determinantal 
point process in M.'^. First, we make the simple observation that the pair correlation 
function specified by ([T]) is related to the Fourier transform a{k] R) of the one- 
dimensional scaled intersection volume evaluated at = r and R = it, namely, 

92{r) = 1 - (36) 

where a{k; R) is given by for d = 1. A natural generalization of this pair correlation 
function in R'^ is to replace the one-dimensional Fourier transform in (136!) with its d- 
dimensional counterpart (j24l) evaluated at = r and divided by (27r)'^, i.e., 

a{r; R) 

However, the value of R in each dimension must be chosen so that such a pair correlation 
function corresponds to a determinantal point process in that dimension. In other words, 
if we take the positive semidefinite function H{r) to be given by the following radial 
function: 

_ ^/W^ _ Vni + d/2)J,/2{rR) 

R must be determined so that the conditions fl2^ and (1281) are satisfied. Noting that 
the expansion of \H{r)\'^ for small r to leading order is given by 

l^"-'l'- 2W(lV./2) +°<'''' P^' 

and using the condition that (72 (0) = [cf. ( l28i) ] yields that 

R = K = 2y^[T{l + d/2)Y^'^. (40) 

Now we must show that a pair correlation function (137|) with R = K satisfies the 
bounds of (1281) and the bounds (l26ll on the spectral function H{k). The square of the 
function H{r), specified by f l38|) . is clearly positive and achieves its maximum value of 
unity at the origin when R = K and tends to zero in the limit r — 00, and hence the 
bounds of (!28|) are satisfied when R = K. Referring to relation (123!) for the Fourier 
transform of the window indicator function, we see that the Fourier transform of (!38ll 
with R = K is simply the indicator function 

H{k) =w{k;K) = e{K -k), (41) 
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which automatically satisfies where w is specified by ( IT^ and K is given by (HO]) . 

In summary, we have demonstrated that there is a determinantal point process in 
M*^ with n-particle densities given by ( l25l) with the kernel 



where i^' is given by ( l40l) . We will call such a determinantal point process a "Fermi- 
sphere" point process because H{k) corresponds to a spherical window indicator function 
in Fourier space and, as we will see, corresponds to the ground state of a fermionic system 
in which the Fermi sphere is completely filled (see Section 3). In particular, the pair 
correlation function of such a point process [33] is given by 

,,(r) = l-2^r(l + d/2f^|-^. (43) 

Moreover, the corresponding structure factor, at unit number density, takes the form 

S{k) = l- a{k;K), (44) 

where a{k; K) is the scaled intersection volume of two (i-dimensional spheres of radius 
K separated by a distance fc, i.e., the function a(r;_R), specified by (I2T1) . with the 
replacements r k and R K. Relation (l44l) is easily obtained by taking the Fourier 
transform of the total correlation function defined by (1291) . where H{r) is given by fH2]) . 
and employing the definition ([7]) for the structure factor. It follows from the properties 
of a{k; K) in (jSj) that the structure factor S{k) obeys the following bounds for all k: 

< S{k) < 1, (45) 

and achieves its maximum value of unity for k > 2K. The corresponding cumulative 
coordination number is given by 

Z{r) =vi{r)-d / -dx. (46) 

./n X 



Note that the first term in ( 1461) is precisely the cumulative coordination number for 
a Poisson point process in and that the second term is strictly negative for any 
r > 0, which is a refiection of the short-range repulsive nature of the point process. We 
will show that for sufficiently large r, the dominant contribution to Z{r) will be the 
Poissonian term w(r). We will see it is the cumulative coordination number Z{r), rather 
than the pair correlation function (contrary to conventional wisdom), that enables one 
to identify and quantify an effective "hard-core" diameter of the determinantal point 
processes. 

It is instructive to examine the asymptotic behaviors of g2{f), S{k), and Z{r) for 
the Fermi-sphere point process. By virtue of the asymptotic properties of the Bessel 
function of arbitrary order, the small-r and large-r forms of the pair correlation function 
( l43l) are respectively given by 
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and 



,, ^ 2T{l + d/2)cos^rK -7r{d+l)/4) . . 



We see that g2{'f) tends to zero quadratically in r in the hmit r — 0, independent 
of the dimension. The coefficient of the quadratic term in (1471) for positive d attains 
it maximum value of 7r^/3 = 3.2898 . . . for d = 1 and monotonically decreases with 
increasing dimension, asymptoting to 27r/e = 2.3114 ... in the limit d oo. Moreover, 
g2{r) tends to unity for large pair distances with a decay rate that is controlled by the 
power law l/r'^^-^ for any d >1. 

The latter result implies that the structure factor S{k) tends to zero linearly in k 
in the limit k 0. Using ( |2Ti) and (jSj), it is easy to verify that as k tends to zero at 
p=l 

S{k) = '^k+ Oik') (k^O), (49) 

where c{d) is a d-dependent positive constant given by fl20l) . We see that at unit density, 
the point process in M*^ for any d is hyperuniform, and the structure factor S{k) has a 
nonanalytic behavior at the origin given by S{k) ~ \k\ {k ^ 0). The coefficient of the 
linear term in (1491) for positive d attains it minimum value of l/(27r) = 0.1591. . . for 
d = 1 and monotonically increases with increasing dimension such that it asymptotes 
to e^/^/(27r) = 0.2624. . . in the limit d oo. In three dimensions, this unusual linear 
nonanalytic behavior of the structure factor at = is a well-known property of bosons 
in their grounds states [lOl [TTl [12] and, more recently, has been shown to characterize 
maximally random jammed sphere packings [T3] . 

In Fig. [3l we compare pair statistics in both real and Fourier space for d = 1 and 
d = 3 at unit density. We see that the amplitudes of the oscillations in g2{r) that are 
apparent for d = 1 are significantly reduced in the corresponding three-dimensional pair 
correlation function. The smallest value of r for which g2{r) attains its maximum value 
of unity, which we denote by tq, is determined by the ffist positive zero of Jci/2{Kr), 
which for sufficiently large d is given by 

d . 1-301687 0.007942 / 1 \ 
Kro = - + 1.4729154^1/^ + ^ + ^ ( rfV^ ) (d^^)- (50) 



Since K increases with increasing d and grows like y'2nd/e for large d, ro grows like the 
square root of d for large d. Similarly, the smallest value of k for which S{k) attains its 
maximum value of unity grows with increasing d and for large d grows like y/d. 

The results immediately above enable us to obtain the following small-r and large-r 
forms of the cumulative coordination number: 



Z{r) = Vi{r) 
and 



dK^ 2 d(d + 3)K^ 4 , 



0) (51) 



_(ci + 2)2 2(d + 2)2(ci + 4)2 
Z{r) =vi{r) -l + 0{r-^) (r ^ oo). (52) 
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Figure 3. Left panel: Comparison of the pair correlation functions for d = 1 and 
d — 3. Right panel: The corresponding structure factors. 



It is straightforward to prove that the asymptotic result ( 152|) is a rigorous lower 
bound on Z{r) for all r using the identity dxJ'^i^^xK) / x = 1/d, the inequality 
dxJli2{xK) /x > dxJ^i^^xK) /x for all positive r, and relation (l46ll . Indeed, since 
Z{r) is always nonnegative, we have the lower bounds 




Z(r)>{ , '^^^ (53) 



where the length scale 



D = £(l±^ (54) 



is the zero of fi(r) — 1. The length scale D can be regarded as an upper bound on the 
effective hard-core diameter, which clearly grows with increasing dimension as vividly 
illustrated in Fig. |U which compares the exact form of Z[r) with the lower bound fl53l) 
for (i = 3 and d = 17. For large d, the effective hard-core diameter is given by the 
asymptotic expression 



D ' ' 



i + !f^ + o(i) 



(d^oo), (55) 



27re 

which is seen to grow like the square root of d. This growth of the effective hard- 
core diameter with dimension is a conceptually important conclusion that runs counter 
to conventional understanding of corresponding fermionic systems, which identifies the 
inverse "Fermi" radius, i.e., = l/[2-y7rr(l + 1/dY^'^] (a decreasing function of d) 
with an effective hard-core diameter. We elaborate on this point in Section 3. 

The hyperuniformity of the point process implies that the number variance cr^(i?) 
for a spherical window of radius R must grow slower than (i.e., the window volume) 
for large R. However, the fact that the large-r behavior of the pair correlation function 
is controlled by the power law means that cr^(i?) must also grow faster than the 




Figure 4. Cumulative coordination number Z(r) as a function of r for d = 3 and 
d = 17 compared to the corresponding lower bounds obtained from the inequality 



surface area of the window oi ^ . In particular, upon substitution of (H3|) into (fT3|l . 
an asymptotic analysis reveals that for large R 

-^(«' = { 2r|(d+iV2rrr+^/2)-/- '°<^>^^''''}^" 

(i? ^ oo), (56) 

where C{d) is a (i-dimensional constant of order unity. We remark that a similar 
asymptotic scaling is expected to hold even when the observation window is non- 
spherical; a discussion of this point has been provided in [M]. We see that the number 
variance scaled by the window surface area, a"^ (R) / R'^~^ , grows like ln(i?), independent 
of the dimension. This unusual number variance growth law in three dimensions has 
also been recently seen in maximally random jammed sphere packings |T3j, which 
can be viewed as prototypical glasses because they are simultaneously perfectly rigid 
mechanically and maximally disordered. Note that the coefficient multiplying ln(i?) in 
( !56l) decays to zero exponentially fast as ^ oo, and, therefore, the surface- area term 
j^d-i increasingly becomes the dominant one in the large-c? limit. This behavior should 
be contrasted with number variance for a Poisson point process, which grows like the 
window volume, i.e., R'^. 

Figure [5] graphically depicts the Fermi-sphere point processes in one and two 
dimensions, which are generated using the algorithm of Hough et. al. [29]. Details 
and applications of this algorithm are reported by us elsewhere [H]. 
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Figure 5. Top panel: A realization of 50 points of a Fermi-sphere point process in a 
linear "box" subject to periodic boundary conditions. Bottom panel: A realization of 
100 points of a Fermi-sphere point process in a square box subject to periodic boundary 
conditions. 



4-2. "Fermi-Shells" Point Processes in 



Here we consider a generalization of the Fermi-sphere point process in which H{k) is 
an indicator function for concentric rings in reciprocal space; we denote the resulting 
determinantal process as the "Fermi-shells" point process. Without loss of generality, 
define 2m radii kp\ where j G {1, 2, . . . , 2m — 1, 2m} and m G N, such that the region of 



(i-dimensional space within the ball -8(0; kp^) is empty, the annulus between kp' and k 
is filled, and so forth. We therefore have m filled concentric rings in reciprocal space. 
Since the conditions for a determinantal point process are fulfilled for any indicator 
function in reciprocal space [35], the pair correlation function will still be given by (J2711 . 
The calculation of H{r) proceeds as follows: 

d m 



,(2) 



.(2i-l) 



k)} 



(57) 



which implies: 



d/2 m 



d/2 



.(2i-l) ' 



d/2 



■Jd/2ik 



{2j-l) 



(58) 



where ^ denotes the Fourier transform to coordinate space. 
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It is important to note that the values of the various kp^ are not independent of 



each other and are constrained by the density (here set to unity) according to: 



E 

j=i 



(2i) 



,(2j-l) 



i2V^Yr{i + d/2). 



(59) 



The filhng of Fermi shells generally introduces a greater level of short-range correlations 
relative to the Fermi-sphere case (see Fig. [6]). 
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Figure 6. Comparison of the pair correlation function of the Fermi-sphere point 
process (PP) to that of the Fermi-shells point process for d = 2. In the latter case, 



kp = k'p'^ = 4 and k'^'^ = 



kl 



{k\ 



4-3. Spin- Polarized Fermionic Gas in 

Note that we do not know of any correspondence of the general Fermi-shells point process 
in M.'^ for any d > 2 to random matrix theory or the zeros of any generalized zeta function 
in number theory. However, we can show that a spin-polarized fermionic gas in M*^ for 
any d > 1 has the same n-particle densities. For simplicity, this comparison will be done 
for the Fermi-sphere case, i.e., we will demonstrate that the pn are given by fl25|) with 
H{r) given by provided that the Fermi radius kp = K = 2v^r(l + d/2Y''^. 

We first recall some general properties of the n-particle density functions p„ for 
a ground-state noninter acting gas of fermions in W^,d > 1. Assume that we have N 
spin-polarized fermions in a volume V with number density p = N/V. The n-particle 
density functions are defined for all n < with respect to the ground state \(f)o) by: 

n 

p„(Xi,...,X„) = (0o| n[^*(^»)^(^»)] l^^o)' (60) 

i=l 

where the operators '?/'*(x), ?/'(x) are the creation and annihilation field operators, 
respectively. Using the momentum representation, the field operators are defined in 
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terms of the particle creation and annihilation operators aj^, a^: 

?/'(x) = ^ 0k(x)ak (61) 

k 

r{^) = J2^< (62) 



k 



1 



1/2 



.k(x)=(-) exp[i(k,x)] (63) 



1 



1/2 



0k(x)=(^-J exp[-i(k,x)], (64) 

where (k, x) = Yl'i=i ^i^i is the Euclidean inner product of two real- valued vectors. The 
field operators for fermions must satisfy the following anticommutation relation: 

{^(x,),V^*(x,)} = 5(x,-x,). (65) 

Wick's Theorem along with the anticommutation relation (I65p immediately allow 
us to write the expectation value in fl60l) as a determinant (equivalently, see Macchi's 
discussion of fermion processes p7]): 



p„(xi,...,xj = ^^(-1)'^ J](0o|^*(x.)^(x.(,))|0o) (66) 

o-G5„ i=l 

= det[(0o|^*(x.)^(x,)|</>o)]i<,,,-<„, (67) 

where Sn denotes the permutation group for n objects. Namely, Wick's Theorem relates 
the n-particle density functions to the n\ ways of pairing the 2n creation and annihilation 
field operators; anticommutation of the operators introduces the factor of (—1)°^ in ( |66|) . 

For the ground-state system, the argument of the determinant in (1671) may be 
evaluated by filling the Fermi sphere up to the Fermi radius kp with the result: 

G(x,x') = (0o|^*(x)^(x')|0o) = l^] r-^/V,/2(A;^|x-x'|). (68) 

Since G is a translationally-invariant function, we will abuse notation slightly and write 
G(x, x') = G(|x — x'l) = G{r; kp), making the parameterization by kp explicit. 

In accordance with our established convention and without loss of generality, let 
p = 1. For a system of N spin-polarized fermions, the Fermi radius kp is given exactly 
by: 

kp = 2v^ {F [(rf/2) + 1]}'^" = K. (69) 

Therefore, we may equivalently write for (1681) : 

G(r; K) = TT-'/yr [(rf/2) + l]r-^/Vrf/2(rir) = H{r). (70) 

The result in (1701) shows that our proposed Fermi-sphere point process corresponds 
exactly to the one generated by a system of noninteracting fermions in d dimensions. 

The connection for d = 1 between this system of noninteracting fermions and the 
CUE of random matrix theory implies that the correlations in the ground state resulting 
from the Pauli exclusion principle are equivalent to those induced on a = 2 Coulomb 
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gas constrained to the unit ring and interacting via a logarithmic pair potential at finite 
temperature T = 1/2 in units such that Boltzmann's constant fc^ = 1. The argument in 
Appendix A shows that this reduction of the probability density to a classical particle 
system with at most two-body interactions is peculiar to the choice of the indicator 
function in reciprocal space; in general, one must include at least three-body interactions 
to describe the system appropriately. 

We also mention the significance of flMl) in defining an effective hard core on 
the system of noninteracting fermions in each dimension. This issue is equivalent 
to assigning an appropriate length scale for the widely-studied (see, e.g., [39]) Fermi 
correlation hole. An argument primarily due to Slater [50] suggests that the correlations 
which result from antisymmetry in the many-body wavefunction extend outward for 
a distance tq ^ kp^. The reasoning behind this choice of length scale relies on the 
introduction of an exchange hole into the system with approximate spherical symmetry 
and a radius determined by the de Broglie wavelength associated with the Fermi radius 
kp- This notion has been quantified for interacting atomic systems by considering, for 
example, the difference of distribution functions derived from Hartree-Fock and Hartree 
wavef unctions [H], but the kp^ scaling for the noninteracting case is still reported in 
modern texts [12] . 

However, we recall from (HOj) that as d increases, the value of the Fermi radius also 
increases; in fact, kp ~ ^/d for sufficiently large d. This behavior implies that kp^ — > 
as d —>■ +00, meaning that the effective hard core diameter would become negligibly 
small for sufficiently large d if it were to scale as kp^. This conclusion contradicts the 
fact that g2{r) pointwise as d ^ +oo for all r G [0, +oo). In other words, the range 
in r over which one finds a small probability of finding two particles in close proximity 
increases with dimension, requiring a different means by which to quantify the extent 
of the effective hard core. 

The definition of Z{r) in suggests instead that we take the value of D in ( 15^ as 
a measure of the effective hard core diameter. We note that since Z{r) is a nonnegative 
monotonically increasing function of r and, therefore, nonzero for some range in r over 
which its lower bound is zero, D must represent an upper bound to the effective hard 
core radius. This representation of the effective hard-core diameter is quantitatively 
well-defined for any dimension due to the inclusion of the x'^~^ factor under the integral 
in ([8]) from the surface area of the d-dimensional ball; it is this factor which appropriately 
rescales g2 such that (l53i) represents a true lower bound on Z{r) in any dimension by the 
argument above. It is also for this reason that neither (72 (r) nor S{k) alone are sufficient 
to define the effective hard-core diameter in a quantitative manner. In contrast to 
Slater's scaling of the effective hard core diameter, D ^ kp in any dimension by the 
definition of the Fermi radius, which is in accordance with the high- dimensional behavior 
of both g2 and kp. 
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5. Nearest-Neighbor Functions 

It is useful to characterize point processes by examining other quantities besides the 
ra-particle correlation functions. One popular descriptor used in one dimension is the 
so-called "gap" distribution function p{z) [3], which characterizes the spacing between 
the points. In the random matrix theory literature, this quantity often has erroneously 
and misleadingly been called the "nearest-neighbor-spacing" distribution because gaps 
to the right of some reference point are considered. However, p{z) makes no distinction 
between gaps to the left or right of some reference point. The quantity p{z)dz gives the 
probability of finding a gap (a line interval empty of points) of length between z and 
z + dz. The function p{z) is called the chord length probability density in the theory of 
random media [361 EZ] • 

In the case of random matrix theory, there exist exact representations of p{z) for 
the spacings of the eigenvalues in the GOE, GUE, and GSE, but they can only be 
determined numerically for general situations. A remarkably accurate approximation 
for the GOE in the infinitely-large matrix limit is the so-called Wigner surmise. The 
Wigner surmise has been generalized to any of the aforementioned ensembles as follows: 



where the parameters A/^ and B/j, which depend on reciprocal temperature /3, are 
obtained from the normalization of both p{z) and its first moment, or the average 
gap size (z). For the GOE, GUE, and GSE, /? = 1, 2, and 4, respectively, and Ai = n/2, 
Bi = 7r/4, A2 = 32/7r2, B2 = 4/7r, and A^ = 262144/(7297r2), B^ = 64/(97r). 

For d > 2, the gap distribution function is strictly not a meaningful descriptor 
of point processes. The natural generalizations of p{z) in higher dimensions are the 
nearest-neighbor functions [3H1 ETj . Nearest-neighbor functions describe the probability 
of finding the nearest point of a point process in at some given distance from a 
reference point in the space. Such statistical quantities are called "void" or "particle" 
nearest-neighbor functions if the the reference point is an arbitrary point of the space 
or an actual point of the point process, respectively. 

5.1. Definitions 

First, we recall the definitions of the void and particle nearest-neighbor probability 
density functions Hv{r) and Hp{r), respectively: 

Hv{r) dr = Probability that a point of the point process lies at 

a distance between r and r + dr from an arbitrary (72) 
point in the space. 



pp{z) = Apz^e-""^ 



(71) 



Hp{r) dr 



Probability that a point of the point process lies at 

a distance between r and r + dr from another point (73) 

of the point process. 
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It is useful to introduce the associated dimensionless "exclusion" probabilities Ey{r) 
and Ep{r) defined as follows: 

Ey{r) = Probability of finding a spherical cavity of radius r 
empty of any points in the point process. 

Ep{r) = Probability of finding a spherical cavity of radius r 

centered at an arbitrary point of the point process (75) 
empty of any other points. 
It follows that the exclusion probabilities are complementary cumulative distribution 
functions associated with the density functions and thus are related to the latter via 

/•r 

Evir) = l- / Hvix)dx (76) 
Jo 

and 

/•r 

Ep{r) = l- / Hp{x) dx. (77) 



Differentiating the exclusion-probability relations with respect to r gives 

»v(r) = (78) 

and 

HAr) = (79) 
The nth moment of Hp{r) is defined as 

POO 

\n= r''Hp{r)dr. (80) 
Jo 

Of particular interest to us is the mean nearest-neighbor distance 

POO 

A = Ai = / rHp[r)dr 



Ep{r)dr. (81) 

It is useful to express the density functions Hy{r) and Hp{r) as a product of two 
functions as follows: 

Hy{r)=ps^{r)Gv{r)Ey{r), (82) 

and 

Hp{r)=ps^{r)Gp{r)Ep{r), (83) 

where si(r) is the surface area of a d-dimensional sphere of radius r given by ([9]). The 
quantities Gy{r) and Gp{r) are called the "conditional" nearest-neighbor functions and 
have the following interpretation: 

psi{r)Gv{r)dr = Given that a spherical cavity of radius r 
centered at an arbitrary point in the space is 
empty of any points of the point process, the 
probability of finding a point in the spherical 
shell of volume Si(r)dr surrounding the 
arbitrary point. 
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psi{r)Gp{r)dr = Given that a spherical cavity of radius r 
centered at a randomly selected point of 
the point process is empty of any other 
points, the probability of finding a point 
in the spherical shell of volume Si(r)dr 
surrounding the randomly selected point. 

Thus, it follows that the exclusion probabilities are also given by 



and 



Ev{r) 



Epir) 



exp 



exp 



-psi{l) 



X 



d-1 



Gp{x)dx 



i5) 



(86) 



^7) 



We remark that knowledge of any one function H, E, or G (either void or particle) is 
sufficient to determine the other two functions via the relations mentioned above. 



5.2. Series Representations 

The nearest-neighbor functions can be expressed as infinite series whose terms are 
integrals over the n-particle density functions [3HI ETj. For example, the void and 
particle exclusion probability functions for a translationally invariant point process are 
respectively given by 

oo k r ^ 

Ev{r) = 1 + E(-l)'lT / 9k{r') H ©(^ " 1^ " r.l)dr, (88) 



k=l -Ji^- j^i 



and 



fc+i 



Epir) = 1 + E(-l)'fr / 9k+i{r'^') J] ©(^ " ^i.)dr.- (89) 

k=l ' ^^"^ j=2 

The corresponding series for Hv{r) and Hp{r) are obtained from the series above using 
(CHD and (CSD. 

In general, an exact evaluation of the aforementioned infinite series are not possible, 
except for simple processes such as the Poisson point process. In the latter instance, 
because p„ = p", both series flHHj) and fl89|) can be summed exactly to give 

Ev{r) = Ep{r) = exp[-pt;i(r)], (90) 

where f i(r) is the volume of a (i-dimensional sphere of radius r given by ( fT4l) . Therefore, 
for a Poisson point process, we have from (ITHI) . (179|) . (!82l) . and ( !83|) that 

Hy{r) = Hp{r) = psi(r) exp[-pt;i(r)], Gyir) = Gp(r) = 1. (91) 

We see that there is no distinction between the void and particle quantities for the 
Poisson distribution, which is generally not the case for correlated point processes. 
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5.3. Rigorous Bounds 

Torquato has given rigorous upper and lower bounds on the so-called canonical n-point 
correlation function iJ„(x™; x^"*"; r'^) (with n = p + q and m < p) for point processes in 
M*^. Since the void and particle exclusion probabilities and nearest-neighbor probability 
density functions are just special cases of then we also have strict bounds on them 
for such models. Let X represent either Ey, Hy, Ep, or Hp and X^''^ represent the kth 
term of the series for these functions. Furthermore, let 

e 

= ^(-1)^X('=) (92) 

k=0 

(93) 

be the partial sum. Then it follows that for any of the exclusion probabilities or nearest- 
neighbor probability density functions, we have the bounds 

X <W\ for £ even 

X > W\ for i odd. (94) 

Application of the aforementioned inequalities yield the first three successive bounds 
on the nonnegative void exclusion probability: 

Ev{r) < 1 (95) 
Ev{r) > 1 - pvi{r) (96) 

2 r2r 

Ev{r) < 1 - pvi{r) + ^Si(l) / x'^-^<*(x; r)g2{x)dx, (97) 

^ Jo 

where v^^{x;r) = Vi{r)a{x]r) is the intersection volume of two d-dimensional spheres 
(cf. flTSl) ). The corresponding first two nontrivial bounds on the nonnegative nearest- 
neighbor probability density function Hy{r) are as follows: 

Hv{r)<psi{r) (98) 

2 r2r 

Hvir) > psi(r) - ^s,{l) / x''-hl^\x;r)g2{x)dx, (99) 
^ Jo 

where s^^^{x; r) = dv^^{x; r)/dr is the surface area of the intersection volume v^^{x; r). 
Bounds on the conditional function Gy(r) follow by combining the bounds above on 
Ev{r) and Hv{r) and definition (|82|) . For example, we obtain the following bounds 

Gy{r) < (100) 

1 - pvi[r) 



and 



I- jf^jrMl) C X'' H'^\x-r)g2{x)dx 

Gv[r) > — — 5 -1^ ; , (101) 

1 — pvi{r) + ySi(1) Jq x'^~^f2"*(x; r)g2{x)dx 



which should only be applied for r such that Gv{r) remains positive. The bounds above 
lead to the following properties of the nearest-neighbor functions at the origin: 

Ev{0) = 1, Hv{0) = 0, Gv{0) = 1. (102) 
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Similarly, the first three successive bounds on the particle exclusion probability are 
given by 

Ep{r) < 1 (103) 
Ep{r) > 1 - Z{r) (104) 
Epir) < 1 - Z{r) 

+^ / / 0(?^ - n2)0(r - ri3)5(3(ri2,ri3,r23)dr2dr3, (105) 

where Z{r) is the cumulative coordination number defined by ([8]). The corresponding 
first two nontrivial bounds on the nonnegative nearest-neighbor probability density 
function Hp{r) are as follows: 

Hp{r) < psi{r)g2{r) (106) 
Hp{r) > psi{r)g2{r) 

-p^ / S{r -ri2)&{r -ri3)g3{ri2,ri3,r23)dr2dr3, (107) 

JRd JRd 

where 6{r) is the radial Dirac delta function. Bounds on the conditional function Gp{r) 
follow by combining the bounds on Ep{r) and Hp{r) and definition ( l83i) . For example, 
we obtain the following bounds: 

and 

^?2(r) - ^4,4,5(r-ri2)e(r-ri3)(73(ri2,ri3,r23)dr2dr3 

Gp(r) > , (109) 

1 - Z{r) + ^ J^a J^d Qir - ri2)0(r - ris)gs{ri2, r^, r23)dr2dr3 

which should only be applied for r such that Gp{r) remains positive. The bounds above 
lead to the following properties of the nearest-neighbor functions at the origin: 

Ep{0) = 1, Hp{0) = 0, Gp{0) = 0. (110) 

We now obtain bounds on the mean nearest-neighbor distance A at unit density 
using the aforementioned upper and lower bounds on Ep{r). Let us define the following 
distances: 

\l= / [1-Z(r)]dr (111) 
Jo 

and 

Xu= / exp[-Z(r)]dr, (112) 



where r^ is the location of the zero of 1 — Z{r). In light of the bounds (I104p and 
it is clear that and Xu bound A from below and above, respectively, i.e., 

Al < A < Xu. (113) 
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5.4- Results for Fermi-Sphere Point Processes 

5.4-1 ■ Exact Determinantal Representations The Fermi-sphere point process is unique 
in that both the particle and void exclusion probabilities may be expressed as 
determinants over N x N matrices, the elements of which are related to overlap integrals 
of the basis functions 0k = (l/VV) exp[i(k, x)] on -8(0; r), a d-dimensional ball of radius 
r centered at the origin (the exact location of the ball's center is irrelevant since the 
point process is translationally invariant). The thermodynamic limit can then be taken 
appropriately. We provide the details of this analysis elsewhere [H] and immediately 
state the results: 

Ey(r) = det [I- M(r)] (114) 
Ep{r) = Ev{r) tr{A[I - M(r)]-^}, (115) 

where I is the N x N identity matrix, and the matrices M{r) and A are defined by: 



M,,(r) = / 0,(x)0,(x)dx (116) 

JB{0;r) 

A, = 0l(O)0.(O)/p, (117) 

where p is the number density. We recall that knowledge of Ev{r) and Ep{r) is sufficient 
to determine all of the remaining nearest-neighbor functions Hy/p{r) and Gy/pir). Note 
that as r — > 0, Mij{r) —>■ for all i and j, which provides the necessary results Ev{0) = 1 
and Ep{0) = ti^A) = J2i l0i(O)P/P = -^(0) = 1- We mention that these determinants 
must be evaluated numerically but essentially exactly for finite N x N matrices, where 
N is chosen to be sufficiently large to capture accurately the behavior of the system in 
the thermodynamic limit. Evidence for this convergence is provided in another paper 

m 



5.4-2. Bounds, Comparison to Exact Results, and Link to Sphere Packings We now 
obtain bounds for nearest-neighbor functions for the Fermi-sphere point process and 
compare them to the corresponding aforementioned exact results. Using the identity 



si{l) x''+\^\x-r)dx = j^ry,{r)]\ (118) 

the leading order term of the small-r expansion fH7|) . and the upper bound 0971) . we 
obtain the weaker upper bound 

dK'^ 

Ev{r) < 1 - pviir) + p^-^^-—^r^[v,(r)]^, (119) 
which is exact through terms of order r^('^+^). Therefore, we also have 

Hvir) = ps.ir) - P' ^^^J^lf r%s,{r)r + 0{r^'^') (120) 

and 

Gy{r) = 1 + pvi{r) + p\,{rf + ^(r^^+i). (121) 
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We see that through order r^°'+^, G'y(r) > 1. This bound will be shown to apply for 
any r. 

Employing the inequality 5'3('^i2, ?^i3, ^^23) < 5'2('^i2)5'2('^i3) [cf. (EI])] in the upper 
bound (11051) on Ev{r) and lower bound (I107p on Hp{r), we find the following 
corresponding weaker bounds: 

Ep{r) < 1 - Z{r) + ^ (122) 

and 

Hp{r) > psi{r)g2{r) - pSi{r)g2{r)Z{r). (123) 

These bounds in conjunction with the analogous evaluations of the bounds (llOSp and 
(1109P on the conditional pair function Gp{r) yields its exact small-r behavior up through 
terms of order r^: 

Gp{r) = -^r' - .f'^+y y + 0(r6). (124) 
^ d + 2 2(c/ + 2)2(c/ + 4) ^ ^ ^ ^ ^ 

Comparing this expansion to (HTll reveals that Gp{r) = (72 (r) up through terms of order 

We can also show that 

Gv{r) > 1 for all r (125) 

and 

Gp{r) > g2{r) for all r. (126) 

These results are obtained from definitions (!86|) and flHTl) and the following upper bounds 
on the exclusion probabilities: 

Ev{f) < exp[— /)t>i(r)] for all r (127) 

and 

Ep{r) < exp[-Z(r)] for all r. (128) 

To prove fll27p and (11280 . we recall the series representations of Ey and Ep in ( IHHl) and 
f l89|) . respectively, which we may rewrite in the following more compact form: 



hoo 



Ep,y{r) = l + Y,{-lfEfjy (129) 



fc=i 



fe+i 



Ef = P- I g,^,{v') n e(r - ry)dr,. (131) 

It is important to note that the series in fll29p converge absolutely for all r G M''", which 
is easily seen from the inequalities 
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4'^ ^ ^ n^2(ri,)e(r - ri,)dr, = (133) 

Equations (11321) and (I133p follow directly from the inequalities in ( l30l) and (J3T1) . 
respectively. It is therefore true that 

+00 +00 r / x-j^ 

|i?y(r)| < 1 + 5^4'^ < 5^ ^^^^ = expK(r)] < +oo (134) 



k=l k=0 
hoo +00 rr^/ 



|i?p(r)| < 1 + 5^ 4'^ < 5^ = exp[Z(r)] < +oo, (135) 

k=l k=0 



and absolute convergence of the series in (11291) implies convergence of those series, which 

(k) 

in turn allows us to conclude that the sequences Epjy — as /c ^ +oo. 

We now wish to compare the series in (11291) with the following series representations 
of the proposed upper limits in ( I127P and ( 1128^ : 

E'y ' {r) = 1 - pvi{r) + (136) 

k=2 



4^")(r) = 1 - Z{r) + ^^jr^- (137) 

fc=2 

Note that the series in (11291) . (I136p . and (11371) agree up to their second terms; the 
alternating series test then implies that the contributions from the remaining terms is 
no greater than the magnitude of the third terms in the series. Equations (I132p and 
(I133p therefore show that the series in (I136P and (I137P in actuality do bound the series 
in (I129p from above, thereby proving the claims. The lower bounds in (I125p and (11260 
immediately follow from monotonicity and positivity of the exponential function along 
with (IMD and (IS 
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1.4 



Figure 7. Left panel: Upper and lower bounds on Ev{r) for d ~ A and p = 1 as 
obtained from (l96|) and (I127|) compared to the corresponding exact evaluation of it. 
Right panel: Upper and lower bounds on Ep{r) for d = A and p = 1 as obtained from 
()104p and (|128|) . respectively, compared to the corresponding exact evaluation of it. 




Figure 8. Upper and lower bounds on Gp{r) for d = 4 and p = 1 as obtained from 
(|104p and (|128p . respectively, compared to the corresponding exact evaluation of it. 

The left panel of Fig. [7] compares the upper bound fll27p and lower bound fl96l) of 
Ey for a Fermi-sphere point process for d = 4 to the corresponding exact evaluation. 
The right panel of the same figure compares the upper and lower bounds on Ep{r) 
for d = 4 as obtained from fll04p and fll28p . respectively, to the corresponding exact 
evaluation of it. The upper bounds on both exclusion probabilities provide reasonable 
estimates of the exact results as compared to the corresponding lower bounds. Figure [8] 
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compares upper and lower bounds on Gp{r) for d = 4 as obtained from fllU4l) and fll28p . 
respectively, to the corresponding exact evaluation of it. Not surprisingly, the bounds 
become better estimators as the dimension increases and therefore can be profitably 
used in high dimensions, where exact evaluations are difficult to obtain. The left panel 
of Fig. [9] shows the upper bound (11271) and lower bound (l96ll for a Fermi-sphere point 
process for d = 17. It is seen that the bounds essentially converge to unity for the 
range < r < 0.8 and are relatively close to one another for 0.8 < r < 1.1. The 
right panel of Fig. [9] depicts the analogous bounds on Ep{r) for d = 17. Note that the 
bounds on Ev{r) behave almost exactly like the bounds on Ep{r) at this value of d. 
This graphically suggests that as d becomes large, the exact expressions for -Ey(r) and 
Ep{r) approach the same step function, which we demonstrate below. 
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Figure 9. Left panel: Upper and lower bounds on Ev{r) for d = 17 and p = 1 as 
obtained from ([96)1 and (I127p . Right panel: Upper and lower bounds on Ep(r) for 
d = 17 and /? = 1 as obtained from p04| and (fT28| . 



We now show that both Ep and Ey tend to the same step function Q{D — r) in 
the \axge-d limit, where D is the estimate of the effective hard-core diameter defined by 
(ISj). We begin by utilizing the following upper bound on Ep{r): 



Ep{r) < 



1, r<D, 
exp[l — fi(r) 



r> D. 



;i38) 



This upper bound is obtained by combining lower bound (153|) on the cumulative 
coordination number Z{r) and upper bound (11281) . and hence is a weaker upper bound 
on Ep{r) than f ll28p . As d tends to infinity, we see that the upper bound (11381) tends 
to the step function 6(D — r), i.e., 

Ep{r)<Q{D-r) {d ^ oo). (139) 



Now we show that Ep is bounded from below by the same step function in this limit. 
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Ep{r)>e{D-r) {d ^ oo). (140) 

To prove this, we note that the the lower bound (11041) on Ep{r) tends to a unit step 
function as d becomes large [cf. Fig. [9] whose discontinuity location cannot exceed the 
zero of 1 — Z(r), which can be estimated to be given by 

\d + 2fT{l + dl2) V'^'^'^ 

dK^ \ ' ^ ' 

This estimate, which bounds the zero from below and becomes increasingly accurate 
as d tends to infinity, is obtained by substituting the leading term of the asymptotic 
expansion (|5T1) into 1 — Z{r) and solving for Dq. For large d, Dq has the asymptotic 
expansion 



, d 



i-!^ + o(i) 



(d^oo), (142) 



27re 

Comparison of this expansion to that of the corresponding one for D [cf. (l55l) ] reveals 
that D bounds Dq from above for sufficiently large d and Dq ^ D in the limit d ^ oo. 
Thus, the lower bound on Ep{r) in this limit is given by (11401) . Combination of this 
lower bound with upper bound (11391) leads to the following high-dimensional behavior: 

Ep{r)^Q{D-r) (rf ^ oo), (143) 

where we recall that D grows like \fd for large d [cf. (155|) ]. Following the analogous 
analysis using the lower bound ( !96l) and upper bound (1127^ on Ev{r), we can show 

Ev{r)^Q{D-r) (rf ^ oo). (144) 

It is noteworthy that this analysis means that the lower bound (l53l) on the 
cumulative coordination number Z{r) becomes exact in the limit d — oo. This in 
turn implies an "effective" pair correlation function (^) ^^^^ tends to the following 
step function as d tends to infinity: 

gl{r) Q{r - D) (d^oo). (145) 

This effective pair correlation function (^Kr) is to be distinguished from the true pair 
correlation function (137|) . which tends to unity for distances beyond the length scale D, 
but remains a quadratic function of r for small r. Because the effective pair correlation 
function is based on the behavior of Z{r), which weights the true 5'2(r) by r'^~^ (due 
to the appearance of Si(r), the surface area of a sphere of radius r), (^) tends to 
a step function in the high-dimensional limit. The fact that the oscillations of g2{r) 
seen in low dimensions (cf. Fig. [3]) effectively vanish in the large-d limit is consistent 
with the decorrelation principle |16], which, roughly speaking, states that unconstrained 
correlations that exist in low dimensions vanish as d tends to infinity, and all higher- 
order correlation functions Qn for n > 3 may be expressed in terms of g2 within some 
small error. 

We have already shown in (1351) that the latter claim is true for any determinantal 
point process. The former claim is seen from the form of (73 iii (I145p . which immediately 
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suggests that asymptotic unconstrained correlations in the Fermi sphere point process 
diminish with respect to increasing dimension d. In other words, g2 flattens at unity 
for sufficiently large r as d becomes large, which implies that long-range correlations 
between any two particles in the system diminish with respect to increasing dimension, 
leaving only the small-r correlations, which extend outward for a greater range in r 
as d increases. This conclusion in conjunction with fl35|) implies that for sufficiently 
large d and for large particle separations, Qn ~ detl = 1, where I is the n x n identity 
matrix. Therefore, all n-particle correlations also diminish for large particle separations 
and large d in accordance with a decorrelation of the system. We remark that the fact 
that g2{r) ^ as r —>■ does not affect the statement of the decorrelation principle for 
the Fermi sphere point process; borrowing the language of quantum mechanics, these 
correlations arise from the constraint of antisymmetry in the many-particle wavefunction 
and therefore must be enforced in any dimension. 

For sufficiently large d this analysis implies that the system reduces to a sphere 
packing with an effective hard-core diameter equal to D. The connection to sphere 
packings implies that the fraction of space (p covered by the spheres at unit number 
density is bounded from above by the following inequality 

0<^^iP/2) = ^. (146) 

Interestingly, Minkowski proved a lower bound on the coverage fraction of the densest 
lattice sphere packings that asymptotically is controlled by 1/2*^ [43j. We remark on the 
significance of this result in Section 6. 



5.4-3. Mean Nearest- Neighbor Distance We now obtain analytical estimates of the 
mean nearest-neighbor distance A at unit density using the general upper and lower 
bounds on A [cf. fillip and (11121) ]. Consider the arithmetic average of (1 11 II) and (I112p : 



(147) 



For low dimensions, the arithmetic average of the upper and lower bounds for the Fermi- 
sphere point process provide reasonable estimates of A, as seen in Tabled] for the first 
four space dimensions, which also includes the corresponding exact results. We see that 
the estimate A captures the nonmonotone dependence of the mean nearest-neighbor 
distance with dimension. Moreover, the table shows that the upper bound becomes the 
dominant contribution to A as d increases and A becomes increasingly accurate as the 
space dimensions increases. 

Note use of the upper bound (11380 enables us to obtain the following weaker but 
analytically solvable upper bound on A: 



A < Ac/ < A 



\(7* 



D 



D 



D 



exp[l — t>i(r)]dr 

r(iMi)e ' 
d 



;i48) 
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Table 1. Comparison of estimates of the mean- nearest neighbor distance A for the 
first four space dimensions of the Fermi-sphere point process at unit density to the 
corresponding "exact" values. 

d Upper Bound Lower Bound Average of Bounds Exact D 

1 0.917808 0.658199 0.788003 0.725728 0.5 

2 0.688071 0.581193 0.634632 0.649823 0.564190 

3 0.670304 0.593975 0.632139 0.654511 0.620350 

4 0.687631 0.625009 0.656320 0.679561 0.670938 



where T{x, a) is the incomplete gamma function. For large d, we have the asymptotic 
expression 



A 



u* 



D 



D 



1 + 



mi) 
d 



+ o 



1 + 



0.5963473622 . . . 



d 



O 



{d oo). 



Moreover, using the lower bound 

"Do 



A > 



[1 - Z{r)]dr 



(149) 



(150) 



where Dq is the zero of 1 — Z{r), and the aforementioned asymptotic analysis of the 
lower bound on Ep{r), yields 

\>Do (d^oo). (151) 

In summary, combination of the bounds 01481) and (11501) and the asymptotic 
expression (11421) . enables us to conclude that the mean nearest-neighbor distance 
approaches the length scale D as d becomes large, i.e., 

(rf^oo), (152) 

which asymptotically grows like the square root of d according to (155!) and, as we 
concluded above, specifies the location of the step- function discontinuity of Ep{r), Ev{r) 
and g2{r) in the large-d limit. Table [1] shows that the length scale D is already an 
accurate estimate of the mean nearest-neighbor distance for d = 4. 

It is noteworthy that the asymptotic mean nearest-neighbor-distance formula (11521) 
is precisely the same as the asymptotic form of the mean nearest-neighbor distance 
of a Poisson point process. The latter for any dimension at unit density is given by 
r(l + 1/d) r(l + (i/2)^/^/-\/7r [37], which in the high-dimensional limit is exactly equal 
to D. The fact that the mean nearest-neighbor distance for the Fermi-sphere point 
process behaves like that of a Poisson point process in the high-dimensional limit is not 
surprising in light of the decorrelation principle |16j . 

Note that the expression for the mean nearest- neighbor distance A(p) for any density 
p can be related to the corresponding quantity A(l) at unit density by the simple scaling 
relation 

A(l) 



Hp) 



p 



l/d ■ 



(153) 
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Figure [TD] shows the mean nearest-neighbor distance as a function of density for various 
dimensions; the cases d = 1 and d = 4 are exact evaluations and the instance (i = 17 is 
obtained from the upper bound prediction fll48p and the scahng relation f llSSp . which 
is expected to be a highly accurate estimate. 
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P 

Figure 10. The mean nearest-neighbor distance A as a function of density p for various 
dimensions. The cases d = 1 and cZ = 4 are exact evaluations and the instance d = 17 
is obtained from the upper bound prediction p48p in conjunction with p53|) . 

5.4-4- Large-r Behavior We conclude this section by making some remarks about the 
conditional nearest-neighbor functions Gv{r) and Gp{r). The fact that the exclusion 
probabilities -Ey(r) and Ep{r) tend to the same step function in the high-dimensional 
limit implies that G'y(r) and Gp{r) have the same large-r behavior as d tends to infinity. 
In fact, our exact evaluations of Gv{r) and Gp{r) for a finite range of r in low dimensions 
indicate that each function becomes linear in r for large r and the ratio Gp{r) / Gv{r) 
tends to unity. Figure [TT] shows our evaluations of both Gv{r) and Gp{r) for the first 
four space dimensions for the range < r < 1.4 (see our companion paper [13] for 
further numerical details). 

In one dimension, one can show that the nearest-neighbor void functions are directly 
related to a radial function that is a solution to a second-order nonlinear differential 
equation [44]. This differential equation can be evaluated exactly for small and large r. 
In particular, using this asymptotic analysis leads to the following large-r behavior for 
(jy(r) for the one-dimensional Fermi-sphere point process: 

2 

Gy{r) = ^r + 0{r-^) (r ^ oo). (154) 

The coefficient 7r^/2 can be compared to the corresponding Wigner surmise prediction 
of IG/tt, which is obtained via the one- dimensional relation that exactly links Ev{r) to 
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Figure 11. Exact calculations of both Gv{r) and Gp{r) for the first four space 
dimensions. 



the gap distribution function p{z) [15], 

Ev{r)= / {z - 2r)p{z)dz, (155) 



J2r 

and Eq. fl7T]) with j3 = 2. Presumably, Gp{r) has the same asymptotic behavior 
for (i = 1, a conclusion supported by our numerical results [H]. Interestingly, for both 
conditional probability functions to have the same linear behavior for large r implies that 
the corresponding exclusion probability functions have the following large-r behavior for 
finite d: 

Ev{r) = Ep{r) exp[-K(rf)r'^+^] (r ^ cx)), (156) 

where n{d) is a positive d-dependent constant. 

Thus, this analysis reveals that the probability of finding a large spherical cavity 
of radius r for a Fermi-sphere point process in dimension d behaves similar to that of a 
Poisson point process but in dimension d + 1. For a Poisson point process, the constant 
K{d) = 7r('^+^)/^/r[((i + 3)/2], [cf. fl90l) and therefore if Ev{r) for a Fermi-sphere point 
process in M°' for large r behaved exactly like that for a Poisson point process in M'^"'"^, 
then Gv{r) would be given by 



Gv(r)=Gp(r) = ^?||3^, + 0(l) (r ^ oo). (157) 

For the first four space dimensions, the coefficients multiplying r in fll57p are given by 
vr = 3.14 . . ., 2, 7r/2 = 1.57 . . ., and 4/3 = 1.33 .. ., respectively. This can be compared 
to the exact coefficients, which for d = 1,2,3 and 4 are given by 7r^/2 = 4.93..., 
2.815 ± 0.048, 1.597 ± 0.043, and 1.297 ± 0.045, respectively. The last three values 
are estimates that we have obtained based upon extrapolations from the evaluations 
of Gy(r) (see Ref. [H] for further numerical details). We see that the probability of 
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finding a large spherical cavity of radius r in M*^ is approximately the same as that for a 
Poisson point process in M'^+^, implying linear growth of Gv{r) or Gp{r) for large r. It 
is not unreasoanble to conclude that this approximation becomes increasingly accurate 
as d increases. Further justification for this remarkable behavior is given in Section 6. 

6. Conclusions and Discussion 

We have obtained and characterized a new class of determinantal point processes in 
W^, the most general of which we call Fermi-shells point processes. The n-particle 
correlation functions for any n and d, which completely characterize the point process, 
are determined analytically. We focused primarily on a special case, the Fermi-sphere 
point process, which in one dimension is identical to the point process that characterizes 
the spacings of the eigenvalues of the GUE (as well as the CUE), the conjectured spacings 
of the nontrivial zeros of the Riemann zeta function, and the positions of spin-polarized 
fermions in their ground state (i.e., completely filling the Fermi "sphere"). We are not 
aware of any correspondence of the general Fermi-shells point process in M.'^ for any 
d > 2 to random matrix theory or the zeros of any generalized zeta function in number 
theory. If our determinantal point processes have connections to random matrices in 
arbitrary space dimension, the latter must be non-Hermitian. For example, Ginibre [16] 
showed that there are two-dimensional determinantal point processes that correspond 
to complex eigenvalues of random non-Hermitian matrices. Thus, it is an open question 
whether the Fermi-shells point process has any correspondence to random matrices. 

We analyzed in great detail properties of pair statistics, including the pair 
correlation function, the structure factor, and cumulative coordination number, as 
a function of spatial dimension d. The point processes for any d are shown to be 
hyperuniform such that the structure factor (or power spectrum) S{k) has a nonanalytic 
behavior at the origin given by S{k) ~ \k\ {k 0). The latter result implies that the 
pair correlation g2{r) tends to unity for large pair distances with a decay rate that is 
controlled by the power law 1/r'^^^. In three dimensions, such a dominant power- law 
decay of g2{r) is a well-known property of bosonic systems in their ground states p^[T2] 
and, more recently, has been shown to characterize maximally random jammed sphere 
packings [T3j. We also graphically displayed one- and two-dimensional realizations of 
the point processes in order to vividly reveal their "repulsive" nature and demonstrated 
that they can be characterized by an effective "hard-core" diameter that grows like the 
square root of d. 

Our study of the nearest-neighbor functions of the Fermi-sphere point process 
resulted in some noteworthy conclusions. For example, we have seen that the probability 
of finding a large spherical cavity of radius r in M'^ is approximately the same as that for a 
Poisson point process in M"^"*"^, implying linear growth of Gy(r) or Gp{r) for large r. This 
is a remarkable result because it represents the first class of nontrivial point processes 
that we are aware of whose conditional nearest-neighbor functions do not asymptote to 
a constant; see Refs. [3H1 113 HE] and [15] for examples of correlated equilibrium and 
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nonequilibrium point processes that possess Gy(r) or Gp{r) with constant asymptotes, 
respectively. It is clear that the unusual property of linear growth of Gy or Gp for large r 
is due to the long-range nature of the repulsive interactions. Since a Fermi-sphere point 
process in is always hyperuniform, i.e., large wavelength density fluctuations vanish, 
the probability of finding a large spherical cavity must be smaller than the corresponding 
quantity for any point process that is not hyperuniform, such as a Poisson point process 
in W^. The probability of finding a large spherical cavity is of course smaller in a Poisson 
point process in dimension d+1 compared to that in dimension d. Moreover, it is easy 
to show that the asymptotic form (I156P for Ev{r) or Ep{r) for large r is always between 
the aforementioned corresponding rigorous upper and lower bounds on the exclusion 
probabilities. 

We also found that the Fermi-sphere point process becomes a sphere packing in the 
high-dimensional limit with an effective hard-core diameter equal to the length scale D 
[cf. ( 15^ ]. Thus, the fraction of space covered by the spheres at unit number density 
is bounded from above by 1/2'^. This coverage fraction has a special significance in the 
study of sphere packings; it arises not only in Minkowski's famous century-old lower 
bound on the density of the densest lattice sphere packings [13] but in lower bounds for 
saturated and disordered packings [16] as well as the highest achievable density in the 
"ghost" random sequential addition packing [8]. It should be noted, however, that there 
is strong evidence that there exist disordered sphere packings with (p not only greater 
than 1/2*^ [l9l [511 [52|, [50] but with densities that exponentially improve on Minkowski's 
lower bound [16], [23] . 

Elsewhere [H] we report results on the extremes of the nearest-neighbor statistics 
as well as Voronoi statistics of the Fermi-sphere point processes in the first four space 
dimensions. In other work, we will quantify clustering and percolation properties of 
Fermi-sphere point processes. 
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Appendix A. On the Presence of Intrinsic n-Particle Interactions (with 
n >3) for General Determinantal Point Processes 

Our purpose here is to show that the n-particle probability density function for an 
arbitrary determinantal point process, even in one dimension, cannot be written as 
a Boltzmann factor of N classical particles interacting through one- and two-body 
potentials at a finite temperature. Although this claim is true for each of the canonical 
ensembles of random-matrix theory, we show via a counterexample in one dimension 
that intrinsic n-particle interactions with n > 3 are generally necessary to describe the 
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system. 

It is relatively straightforward to express the n-particle probability density function 
for the CUE as a Boltzmann factor of a classical system of pairwise-interacting particles. 
In analogy with the formalism of the one-dimensional Fermi-sphere point process 
introduced in Section 4.1, we fill the Fermi line at a constant density. Thus, may 
write |3j: 

det[exp(inx„)]„,^ = JJ [exp(ix„) - exp(ia;^)]; (A.l) 



n<m 



therefore, 



|det[exp(inxm)]„,„,|^ = exp 



-2 In I exp(ia;„) — exp(ix^ 



exp 



-2 ^ In I sin(a;„ - a;m)/2| 



(A.2) 



(A.3) 



For small eigenvalue separations Aa;„m = 
form: 

I det[exp(ina;m)]n,m|^ ~ exp 



— the result in flA.3P reduces to the 



-2 J]ln|(Ax„^/2)| 



n<m 



(A.4) 



and the probability distribution of the CUE eigenvalues can indeed be written as a 
pairwise interacting potential, which in this example is logarithmic for small particle 
separations. It is important to note that this reformulation of the probability density is 
very peculiar to the method of filling the Fermi sphere. Suppose instead that we decide 
to fill only the states n = 0,2,3, .. .; i.e., the state n = 1 is skipped. We may then 
express the three-body form of the probability density as: 



|det[exp(ina;„)]„,. 
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X sm 



(2vr; 
x-y 



-[3 + 2 cos(a; — y) + 2 cos{x — z) + 2 cos(y — z)] 



sm 



X — z 



sm 



y- z 



(A.5) 



where we have used x,y, and z to represent Xo,X2, and x^; the subscripts denote the 
state n of each particle. The last three factors have the form of the pair interaction in 
OA. 31) . To check whether the pre-factor containing cosines can be written in the same 
form, we write: 

V{x, y,z) = — ln[3 + 2 cos(x — y) + 2 cos{x — z) +2 cos{y — z)]. (A. 6) 
Assume that V{x, y, z) = v{x, y) + v{y, z) + v{z, x). It must then be true that: 



v{x,y) 



1 
2 

ln(3) 



V{x,y,y) - -V{y,y,y) 



3 + + Acos^x - y)], 



(A.7) 
(A.8) 



but by substituting this expression into ([Qj) . we see that we do not recover the 
original functional form. Therefore, we have shown by contradiction that for general 
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determinantal point processes, the interaction potential must contain at least intrinsic 
three-body terms. 
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